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Oscillation frequency of spin torque oscillator with a perpendicularly magnetized free layer and 
an in-plane magnetized pinned layer is theoretically investigated by taking into account the field-like 
torque. It is shown that the field-like torque plays an important role in finding the balance between 
the energy supplied by the spin torque and the dissipation due to the damping, which results in 
a steady precession. The validity of the developed theory is confirmed by performing numerical 
simulations based on the Landau-Lifshitz-Gilbert equation. 


Spin torque oscillator (STO) has attracted much at¬ 
tention as a future nanocommunication device because 
it can produce a large emission power (> 1 /iW), a high 
quality factor (> 10 ^), a high oscillation frequency (> 1 
GHz), a wide frequency tunability (> 3 GHz), and a nar¬ 
row linewidth (< 10^ kHz) [H-Q. In particular, STO with 
a perpendicularly magnetized free layer and an in-plane 
magnetized pinned layer has been developed after the dis¬ 
covery of an enhancement of perpendicular anisotropy of 
GoFeB free layer by attaching MgO capping layer [lOl - 
[H. In the following, we focus on this type of STO. We 
have investigated the oscillation properties of this STO 
both experimentally iEl and theoretically [Ulil. An 
important conclusion derived in these studies was that 
field-like torque is necessary to excite the self-oscillation 
in the absence of an external field, nevertheless the field¬ 
like torque is typically one to two orders of magnitude 
smaller than the spin torque [l 6 l - [l^ . We showed this 
conclusion by performing numerical simulations based on 
the Landau-Lifshitz-Gilbert (LLG) equation [l5l |. 

This paper theoretically proves the reason why the 
field-like torque is necessary to excite the oscillation by 
using the energy balance equation [l9l - lT7j| . An effective 
energy including the effect of the field-like torque is in¬ 
troduced. It is shown that introducing field-like torque 
is crucial in finding the energy balance between the spin 
torque and the damping, and as a result to stabilize a 
steady precession. A good agreement with the LLG sim¬ 
ulation on the current dependence of the oscillation fre¬ 
quency shows the validity of the presented theory. 

The system under consideration is schematically shown 
in Fig. [l](a). The unit vectors pointing in the magnetiza¬ 
tion directions of the free and pinned layers are denoted 
as m and p, respectively. The z-axis is normal to the 
film-plane, whereas the a;-axis is parallel to the pinned 
layer magnetization. The current I is positive when elec¬ 
trons flow from the free layer to the pinned layer. The 
LLG equation of the free layer magnetization m is 


dm 

dt 


ym X H — yiLgin x (p x m) 


7/3iLsm X p -b am x 


dm 
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FIG. 1: (a) Schematic view of the system, (b) Schematic 

views of the contour plot of the effective energy map (dotted), 
Eq. a, and precession trajectory in a steady state with I = 
1.6 mA (solid). 


where 7 is the gyromagnetic ratio. Since the external 
field is assumed to be zero throughout this paper, the 
magnetic field H = (Hk — 47 rM)TOze^ consists of the per¬ 
pendicular anisotropy field only, where Hk and 47 rM are 
the crystalline and shape anisotropy fields, respectively. 
Since we are interested in the perpendicularly magnetized 
free layer, Hk should be larger than AttM. The second 
and third terms on the right-hand-side of Eq. o are the 
spin torque and field-like torque, respectively. The spin 
torque strength, Hg = hr]I/[2e{l -b Am • p)MV], includes 
the saturation magnetization M and volume V of the 
free layer. The spin polarization of the current and the 
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dependence of the spin torque strength on the relative 
angle of the magnetizations are characterized in respec¬ 
tive by rj and A H- According to Ref. , (3 should 
be negative to stabilize the self-oscillation. The values 
of the parameters used in the following calculations are 
M = 1448 emu/c.c., Hk = 20.0 kOe, V = 7rx60x60x2 
nm^, rj = 0.54, \ = , (3 = —0.2, 7 = 1.732 x 10^ 

rad/(Oe-s), and a = 0.005, respectively SEl- The crit¬ 
ical current of the magnetization dynamics for /3 = 0 is 
Ir = [3aeMV/{hriX)\{HY^ — AttM) ~ 1.2 mA, where Ref. 

shows that the effect of (3 on the critical current is 
negligible. When the current magnitude is below the crit¬ 
ical current, the magnetization is stabilized at = 1 . 

In the oscillation state, the energy supplied by the spin 
torque balances the dissipation due to the damping. Usu¬ 
ally, the energy is the magnetic energy density defined as 
E = -M/dm • H [111, which includes the perpendic¬ 
ular anisotropy energy only, —M{Hk — 4:TrM)m1/2, in 
the present model. The first term on the right-hand-side 
of Eq. (IT|) can be expressed as — 7 m x [—dE/d{Min)]. 
However, Eq. (HD indicates that an effective energy den¬ 
sity. 


EeS = 


M{Hk - 47rM) 2 

- 


[3 hr] I 
2e\V 


log (1 -b Am • p), 
( 2 ) 


should be introduced because the first and third terms 
on the right-hand-side of Eq. (HD can be summarized as 
— 7 m X [— i9ifeff/9(Mm)]. Here, we introduce an effective 
magnetic field 'H = —9ifefi/9(Mm) = (/ 3 fi? 7 //[ 2 e(l -b 
Ama,)AfV’], 0, (Mk — 47 rM)mz). Dotted line in Eig. [f] (b) 
schematically shows the contour plot of the effective en¬ 
ergy density Aeff projected to the xy-plane, where the 
constant energy curves slightly shift along the x-axis be¬ 
cause the second term in Eq. (HD breaks the axial sym¬ 
metry of E. Solid line in Eig. HI (b) shows the preces¬ 
sion trajectory of the magnetization in a steady state 
with / = 1.6 mA obtained from the LEG equation. As 
shown, the magnetization steadily precesses practically 
on a constant energy curve of Aeff. Under a given cur¬ 
rent /, the effective energy density Ueff determining the 
constant energy curve of the sta ble p recession is obtained 
by the energy balance equation 


aj!^a{Ees) — ^s{Ees) = 0. (3) 


In this equation, and which are proportional to 
the dissipation due to the damping and energy supplied 
by the spin torque during a precession on the constant 
energy curve, are defined as [l 3 , [ 2 ^ - 1 ^ 


= 7 ^ 


dt 


-{m- ny 


( 4 ) 


f dtHs [p • — (m • p) (m • "H) — ap • (m X 'H)] . 




( 5 ) 

The oscillation frequency on the constant energy curve 


(a) 



(b) 



FIG. 2: Dependences of ^a, —a^a, and their difference 
^a — normalized by ^{Hk — inM) on (0 < rriz < 1) 
for (a) /3 = 0, and (b) /3 = —0.2, where I = 1.6 mA. 

determined by Eq. m is given by 

f = 1/ dt. (6) 

Since we are interested in zero-field oscillation, and from 
the fact that the cross section of STO in experiment Q 
is circle, we neglect external field Hext or with in-plane 
anisotropy field However, the above 

formula can be expanded to system with such effects 
by adding these fields to H and terms —MHext ■ m — 
to the effective energy. 

In the absence of the field-like torque {(3 = 0), i.e., 
Ees = E, there is one-to-one correspondence between the 
energy density E and rriz- Because an experimentally 
measurable quantity is the magnetoresistance propor¬ 
tional to (i?AP — i?p)max[m ■ p] oc maxjma,] = 
it is suitable to calculate Eq. (|3l) as a function of TOz, in¬ 
stead of E, where i?p(AP) is the resistance of STO in the 
(anti)parallel alignment of the magnetizations. Figure H] 
(a) shows dependences of and their differ¬ 
ence .y£a. — oi^a on rriz (0 < < 1) for j3 = 0, where 

and are normalized by 7 (i 7 K — 47rM). The current is 
set as I = 1.6 mA (> 1^). We also show and 

their difference .y£s — aJia for /3 = —0.2 in Fig. (b), 
where uix is set as rrix = — ^/l — Wz- Because —a.Ma 
is proportional to the dissipation due to the damping, 
is always —a.y£a “E 0. The implications of Figs. 
H] (a) and (b) are as follows. In Fig. H] (a), is 

always positive. This means that the energy supplied by 
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current (mA) 

FIG. 3: Current dependences of peak frequency of \'mx{f)\ 
obtained from Eq. ([1} (red circle), and the oscillation fre¬ 
quency estimated by using (|6]) (solid line). 

the spin torque is always larger than the dissipation due 
to the damping, and thus, the net energy absorbed in 
the free layer is positive. Then, starting from the initial 
equilibrium state (ruz = 1), the free layer magnetization 
moves to the in-plane niz = 0, as shown in Ref. [l^. On 
the other hand, in Fig. [2](b), — is positive from 

rUz = 1 to a certain mi,, whereas it is negative from m'^ to 
rUz = 0 (m^ ~ 0.4 in the case of Fig. [l](b)). This means 
that, starting from rUz = 1, the magnetization can move 
to a point m'^ because the net energy absorbed by the 
free layer is positive, which drives the magnetization dy¬ 
namics. However, the magnetization cannot move to the 
film plane {rUz = 0) because the dissipation overcomes 
the energy supplied by the spin torque from rUz = m(, to 


rUz = 0. Then, a stable and large amplitude precession 
is realized on a constant energy curve. 

We confirm the accuracy of the above formula by com¬ 
paring the oscillation frequency estimated by Eq. ® 
with the numerical solution of the LEG equation, Eq. Cl¬ 
in Fig. O we summarize the peak frequency of |ma;(/)| 
for / = 1.2 — 2.0 mA (solid line), where mx{f) is the 
Fourier transformation of rrix (t) . We also show the oscil¬ 
lation frequency estimated from Eq. (jb)) by the dots. A 
quantitatively good agreement is obtained, guaranteeing 
the validity of Eq. ([51. 

In conclusion, we developed a theoretical formula to 
evaluate the zero-field oscillation frequency of STO in 
the presence of the field-like torque. Our approach was 
based on the energy balance equation between the energy 
supplied by the spin torque and the dissipation due to the 
damping. An effective energy density was introduced to 
take into account the effect of the field-like torque. We 
discussed that introducing field-like torque is necessary 
to find the energy balance between the spin torque and 
the damping, which as a result stabilizes a steady preces¬ 
sion. The validity of the developed theory was confirmed 
by performing the numerical simulation, showing a good 
agreement with the present theory. 
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